Exact calculation of the magnetocaloric effect in the spin-1/2 XX Z chain 
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We calculate the entropy and cooling rate of the antiferromagnetic spin-1/2 XX Z chain under 
an adiabatic demagnetization process using the quantum transfer-matrix technique and non-linear 
integral equations. The limiting case of the Ising chain (corresponding to infinitely large anisotropy) 
is presented for comparison. Our exact results for the Heisenberg chain are used as a crosscheck 
for the numerical exact diagonalization as well as Quantum Monte Carlo simulations and allow us 
to benchmark the numerical methods. Close to field-induced quantum phase transitions we observe 
a large magnetocaloric effect. Furthermore, we verify universal low-temperature power laws in the 
cooling rate and entropy, in particular linear scaling of entropy with temperature T in the gapless 
Luttinger-liquid state and scaling as vT at field-induced transitions to gapped phases. 

PACS numbers: 75.10.Pq; 75.30.Sg; 02.70.-c 
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I. INTRODUCTION 

The magnetocaloric effect (MCE) in general addresses 
the change of temperature of magnetic systems under 
the variation of an external magnetic field. The MCE has 
been known since the end of the 19th century^ and it has 
attracted renewed interest recently because of potential 
room-temperature cooling applications (see Refs. 2.3 for 
recent reviews). On the other hand, adiabatic demag- 
netization is a standard low-temperature method: de- 
magnetization of paramagnetic salts was the first method 
to reach temperatures below IK— whereas demagnetiza- 
tion of nuclear spins has reached record low temperatures 
down to lOOpK^ and is still the cooling method of choice 
in the /xK-range.— The cooling rate at the adiabatic de- 
magnetization (§jj) s for an ideal paramagnet (i.e., a 
system of non-interacting magnetic moments) is equal to 
T/H, which means linear monotonic dependence of tem- 
perature on the magnetic field magnitude. The latter is 
a direct consequence of the fact that for any paramag- 
netic system the entropy depends only on the ratio H/T, 
so for any isentrope one gets H/T — const. However, 
the matter could undergo crucial changes for systems of 
interacting spins. For instance, in ferromagnets near the 
Curie point one can observe a substantial enhancement 
of the effect 4 

As has been shown in early investigations, quantum an- 
tiferromagnets are more efficient low-temperature mag- 
netic coolers than ferromagnets i 9 ' 10 ! 11 This fact is con- 
nected with the behavior of the entropy of antiferromag- 
nets. The entropy of any antiferromagnet at low tem- 
peratures displays (at least) one maximum as a function 
of magnetic field, which usually, according to the third 
law of thermodynamics, falls to zero at T — > 0.— The 
Ising model is anomalous in this respect, because of non- 
vanishing zero-temperature entropy at the critical mag- 
netic field H c = qj, where q is the coordination number 
of the lattice and J is the coupling constant. Indeed, un- 



derstanding the influence of quantum fluctuations seems 
to have been an important motivation for the numerical 
work Ref. |T(| The first exact result concerning mag- 
netocaloric properties of the spin-1/2, Ising-like XX Z 
chain has been obtained in Ref. [ill , where the isentropes 
in the (H, T) plane have been presented. The main fea- 
ture of the isentropes of the Ising-like XX Z chain is the 
appearance of two minima. To the best of our knowl- 
edge, the isentropes of the isotropic spin-1/2 Heisenberg 
chain have so far been investigated only numerically 12 
(for numerical studies of the magnetocaloric effect in fer- 
rimagnetic spin chains and higher-spin Heisenberg chains 
see Refs. [l3. 14). Recent measurements^ of the adiabatic 
cooling rate in the spin-1/2 Heisenberg chain compound 
[Cu(/x-C204)(4-aminopyridine) 2 (H 2 0)] n render the mag- 
netocaloric effect of the spin-1/2 Heisenberg chain a topic 
of current interest. 

More generally, the MCE is particularly large in the 
vicinity of quantum critical points (QCPs). The MCE is 
closely related to the generalized Griineisen ratios 
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Here r is the control parameter governing the quantum 
phase transition. In the case of the MCE r is the external 
magnetic field H . Using basic thermodynamic relations, 3 
the generalized Griineisen ratio T h can be related to the 
adiabatic cooling rate (dT/dH), 
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Thus, the magnetic cooling rate is an important quantity 
for the characterization of QCPs, i.e., quantum phase 
transitions between different magnetic structures under 
tuning the magnetic field at T = 0. 

In passing we mention that an analysis of classical 
spin models^ demonstrated that the MCE can be en- 



2 



hanced by geometric frustration. Indeed, adiabatic de- 
magnetization experiments on the frustrated spin-7/2 
pyrochlore-type magnet Gd2Ti 2 07 have shown substan- 
tial drops in temperature in the vicinity of the satura- 
tion field^ Enhanced cooling performance is also theo- 
retically predicted in one-dimensional quantum antifer- 
romagnets such as the J\ — J2-chain and the sawtooth 
chains the diamond chain ) 20 i 21 i 22 ' 23 i 24 as well as in two 
dimensions . 25 i 26 

The one-dimensional spin-1/2 Heisenberg model is fa- 
mous for its integrability. The conventional Bethe ansatz 
technique allows one to obtain all eigenvalues and eigen- 
vectors of the corresponding Hamiltonian, though in non- 
explicit form. There are several sophisticated methods to 
describe thermodynamics of one-dimensional integrable 
models, like thermodynamic Bethe ansatz (TBA), etc. 
(see for example Refs. 27.28j). For pragmatical reasons, 
the most suitable technique is the quantum transfer- 
matrix (QTM) method leading to only two non-linear in- 
tegral equations (NLIEHor the free energy of the Heisen- 
berg chain, see Refs. I29ll30ll3ll and references therein. 
A numerical solution of these NLIEs has been used in 
Ref. |2^ to obtain the magnetic susceptibility and spe- 
cific heat of the isotropic spin-1/2 Heisenberg chain in 
an external magnetic field. However, as far as we are 
aware, there are no exact results for the free energy (or 
equivalently the entropy) and a mixed derivative of the 
free energy which corresponds to the magnetic cooling 
rate in the literature for the isotropic spin-1/2 Heisen- 
berg chain. Filling these gaps is one of the purposes of 
the present paper. 

In this paper we present calculations of quantities re- 
lated to the MCE for the one-dimensional spin-1/2 XXZ 
Heisenberg chain within the QTM and NLIE method. 
Particularly, isentropes for different values of the ex- 
change anisotropy parameter A are obtained. The cool- 
ing rate is computed as a function of external magnetic 
field for various fixed values of temperature and different 
values of the exchange anisotropy. The limiting case of 
infinitely large anisotropy which is just the Ising model 
is considered as well. In that case the expressions for 
all relevant thermodynamic quantities can be obtained 
in closed form. These calculations are supplemented by 
exact diagonalization (ED) and Quantum Monte Carlo 
(QMC) calculations demonstrating full agreement be- 
tween results obtained from the exact solution and nu- 
merical calculations. 



II. ISENTROPES AND COOLING RATE FOR 
THE SPIN-1/2 XX Z HEISENBERG CHAIN 



We will be interested both in the entropy and the asso- 
ciated isentropes as well as the temperature derivatives 
thereof, i.e., the adiabatic cooling rate (j^jj) s - Using 
standard thermodynamic relations one can express the 
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FIG. 1: (Color online) Zero-temperature phase diagram of 
the spin-1/2 XXZ chain in a magnetic field. 



latter as follows^ 
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where Ch is the specific heat at constant magnetic field, 
and M is the magnetization of the system. After normal- 
ization with a factor 1/T, the magnetic cooling rate can 
be identified with the generalized Griineisen ratio Eq. (|2|) 
which we will use in the following. 

We will be specifically interested in the spin-1/2 XXZ 
Heisenberg chain whose Hamiltonian is given by 

N 
n=l 

N 

~Hj2 S n- ( 4 ) 

n=l 

Here N is the total number of sites, S" are spin-1/2 op- 
erators acting at site n, J is the exchange constant, A an 
exchange asymmetry, and H an external magnetic field. 
For the finite-size computations we will assume periodic 
boundary conditions, i.e., S^ +1 = Sf. Note that the 
properties of the XXZ model (U|) are symmetric under 
H -> ~H and M = (S*) -> —M. We will therefore 
concentrate on H > in the following. 

It is useful to recall the zero-temperature phase dia- 
gram of the model Eq. (0}, see Fig. [[]32ii&2i For A > 1 
and small magnetic fields H < H c i there is long-range 
antiferromagnetic order along the z-direction in spin 
space. This state exhibits a gap in the excitation spec- 
trum whose value at H = corresponds to H c \. The 
value of the gap has been computed exactly via the 
Bethe ansatzi 35 i 36 In particular, it is exponentially small 
close to the Heisenberg point at A — 1 (marked by the 
dot in Fig. [T]), which is characteristic for a Kosterlitz- 
Thouless transition , i 38 For A < 1, one finds a gapless 
Luttinger-liquid stated This Luttinger-liquid state ex- 
ists for |A| < 1, H < H c2 and A > 1, H cl < H < H c2 , 
with H C 2 — J (1 + A). Finally, for H > H c2 , the ground 
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state is the ferromagnetically polarized state along the ular the quantum phase transitions at H c \ and H c i will 

^-direction which exhibits again a gap. We anticipate be reflected by the magnetocaloric properties at finite 

that these different zero-temperature regions, in partic- temperature. 

I 



A. Non-linear integral equations 

In this section we recall the approach to the thermodynamics using the QTM technique and NLIEs i 29,31 The 
equations look different for |A| < 1 and |A| > 1. 

First we consider the case |A| < 1. One can represent the free energy of the system per lattice site in the following 
form: 

oo 

f ln[(l + q(x))(l + fl(x))] 
f = e °- T J ^shf dx (5) 

— oo 

(the constant eo is irrelevant in the present context). The auxiliary functions a(x) and a(x) are found from the 
following system of integral equations 

i i s J sin 7 7r irH , . _.. . 

lna(x) = — r2^~ cosh ££ + 2T(n--y) + K *H l + a )( x ) - k + + a)(x) , (6a) 

lna(.) = - J " 2 7 cos l ^ - ^g^y + k * ln(l + a){x) - k~ * ln(l + «)(») ■ (6b) 

Here 7 = arccos A, the symbol * denotes convolution / * g(x) = J"_° f(x — y)g(y)dy and the function k(x) is defined 
by 
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1 7 sinh - ~ 2 fc 
< X ) = TT / V7 7 ' — ^e lkx dk, K ± (x) = K(x±2i). (7) 



2coshfcsinh ( 2 - 1 ) fc 



These equations are valid for < A < 1. Results for negative A can be obtained by changing the sign of the coupling 
J. 

For the case A > 1 the free energy has the form 

tt/2 

/ = e -|c*ln[(l + b)(l + b)] (0) = eo-^ / c(x) In [(1 + b(x))(l + b(x))] dx (8) 

-7T/2 

with / *g(x) = i dy f( x ~ v)9{y) denoting a convolution with modified integration limits and 

00 1 

fc— — 00 

with A = cosh 77. The auxiliary functions are solutions of the following integral equations 

lnb(x) = - JS ™^ c{x) + ^ + k*ln (1 + b(x)) - k+ * In (l + b(x)) (10a) 
In 6(a) = - J c{x) - + A; *ln (l + b») - AT * In (1 + b(a?)) (10b) 

with integration kernels 

fc W= E 9 6 w ^ e ' 2fa ' k ± (x)=k(x±i V *). (11) 
^ z -^ 2cosh(?7K) 

Results for A < — 1 can again be obtained by changing the sign of the coupling J. 
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The equations for the isotropic case A = 1 can be obtained by taking the limit 7^0 from the case |A| < 1 in ([5]) 
and ([5]) or from the case |A| > 1 by changing the spectral parameter x — xrj and taking the limit r\ — > in ([8]) and 

Having all these exact expressions one can obtain any thermodynamic quantity of interest by iteration of the NLIE 
([6]) (or (fT0|) ) and numerical integration of the expression for the free energy ([5]) (or ([8])). Derivatives of the free energy 
with respect to T and H can also be calculated. One can avoid numerical differentiation by solving the associated 
integral equations for the differentiated auxiliary functions, e.g. du lna(a;). Note that derivatives of lno(a;) and lna(x) 
are treated as independent functions in these equations. As an example we will give the equations for the calculation 
of the magnetization per spin M and the derivative of M with respect to the temperature in the regime |A| < 1 



M = T 



4 cosh Tf 



a(x) 
1 + a(x) 



[dn lna(a;) 



a(x) 
1 + a(x) 



[Oh lna(x)] I dx . 



(12) 



The derivatives dn lna(x) and Oh lna(a;) satisfy linear integral equations in which the auxiliary functions a(x) and 
a(x) enter as external functions 



<9jflna(a;) = 
<9fflna(a;) = 



7F7~ T + K * T~r~ $ H lna l ( X )- K * TT= [ q h In a] (x + 2i) , 

21 [tt — 7) 1 + a 1 + a 

- cyrl* \ + K * 7—= [9h In a] (x) - k * [d H In a] (x - 2i) 

21 (7r — 7) 1 + a 1 + a 



(13a) 
(13b) 



To obtain the cooling rate Th we have to determine dAl/dT. In order to achieve this we have to differentiate (|12p with 
respect to T. However it has turned out that in the framework of NLIEs the resulting equations in general behave 
numerically better if the derivatives are taken with respect to the inverse temperature (3 = 1/T (we set ks = !)• 
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a(x) 



[dp In a(x)] [Oh In a(x)] 
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a(x) 
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[dpd H In a(x)] 
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Here four new functions dp\na(x), dpdu fna(x) and their a counterparts occur. The corresponding linear integral 
equations read 



dp lna(x) 
dp lna(x) 
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J sin 7 7r nH 

2 7 cosh if ~ 2 (tt - 7) 
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(15b) 



and 
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[dp In a] [<9 ff In a] + [d^dff In a] 
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Note that the functions dp lna(a;) and dp lna(cc) already allow the calculation of the entropy per spin 

r^r, OO 

ln[(l + a(z))(l + a(z))] , If If a(x) r _ . . , a(x) 



J 4cosh^ T y 4cosh^U + « 

— oo — OO 



[dp In a(x)} + - - [dp In a(x)] J dx . (17) 



For the calculation of the cooling rate Th the specific heat is needed in addition to the mixed derivative dM/dT. It 
can be obtained by differentiating (|17p with respect to /3 and dividing by — T. In the resulting equation another pair 
of functions <9^1na(:r), dp\na(x) occurs, where the corresponding integral equations are derived by differentiating 
([T5|) again with respect to (3 in close analogy to (fl"6)) . 



B. Comparison with numerical results 

First, we present a comparison with numerical meth- 
ods, namely exact diagonalization (ED) and Quantum 
Monte Carlo (QMC). On the one hand, this comparison 
will serve as a cross-check of our results. On the other 
hand, we can use the exact results to assess the perfor- 
mance of the numerical methods. 

The quantities appearing on the r.h.s. of Eq. @ can 
be expressed as follows: 

C H = ^ ((H 2 ) (H) 2 ) , (18) 

{w) = If WH)-{M){H)) , (19) 

where (■) is the expectation value at a fixed tempera- 
ture T and magnetic field H. Here we have chosen a 
normalization per spin which drops out when taking the 
ratio in Eq. ©. Eq. (HHJ) is well known and Eq. (JTSJ) is 
valid for any Hamiltonian conserving magnetization, i.e., 
[M,H] = 0. 

One can write down spectral representations for the 
correlation functions in Eqs. (fT5|) and (fTfJ)) which can be 
evaluated by ED. These correlation functions can also be 
evaluated with QMC. The QMC simulations to be re- 
ported below have been carried out with the ALPS 40 i 41 
directed loop application 4 ^ in the stochastic series expan- 
sion framework^ The specific heat Eq. (fT5|) is measured 
using an improved estimator which involves the fluctua- 
tions of the expansion order 4i The correlation function 
Eq. (fT9|) can be measured in a similar manner. Note that 
it is crucial to choose an appropriate pseudo random- 
number generator in order to obtain correct results. We 
have used the "Mersenne Twister" In our QMC simu- 
lations we have performed 4-10 5 thermalization steps and 
then collected data during a number of sweeps ranging 
between 2 • 10 7 and 3.3 • 10 10 . 

Fig. [5] shows a comparison of Th between the exact 
solution for N = oo, ED for N = 20, and QMC for 
N = 100 at three temperatures which have been cho- 
sen to correspond to the experiments of Ref. [IH The 
ED curves exhibit some wiggles at small magnetic fields 
and temperatures which reflect the fact that the system 
contains only N — 20 sites. In this regime, the QMC 
results for N — 100 are indistinguishable from the ex- 
act solution for N = oo on the scale of Fig. [5] On the 



other hand, the QMC results are subject to big error 
bars at high fields and low temperatures despite the fact 
that we have invested substantial amounts of CPU time 
into these data points. To some extent, this is related 
to performance problems of the algorithm in a magnetic 
field 4^ However, the main reason is that Tjj is given by 
the ratio of two quantities (see Eq. @) which are both 
exponentially small in T for H > H C 2 = 2 J. Indeed, 
even with a large number of sweeps, it is difficult to de- 
termine the ratio of two very small quantities accurately 
by QMC. Conversely, the existence of a gap improves 
finite-size convergence such that ED works particularly 
well in the high-field regime. 

The overall good agreement between all three methods 
serves as a consistency check for each of them. Further- 
more, we see that we can get extremely accurate numer- 
ical results by combining QMC for H < H C 2 and ED for 
H > H c2 (H c2 = 2 J in the present case). 



C. Effects of anisotropy A 

Next, we discuss the effects of exchange anisotropy on 
the cooling rate. Fig. [3] shows results for three selected 
values A = 0.5, 1, and 2. Recall from Fig. Q] that, at 
A = 0.5, one crosses just one quantum phase transition 
at H C 2 = 1.5 J, the curves for A = 1 in Fig. [3] start from 
the special Heisenberg point at H = and cross the tran- 
sition to saturation at H C 2 = 2 J, and the curves at A = 2 
cut two quantum phase transitions, namely H c \ « 0.39 J 
and H C 2 = 3 J. For T = 0.05 J, the quantum phase 
transitions at H C 2 and for A = 2 at H c i are signaled 
by sign changes of the cooling rate Tjj from negative to 
positive values upon increasing field, see Fig.OJa). Note 
that these zeros shift away from the position of the zero- 
temperature phase transitions with increasing tempera- 
ture, as is evident in particular if one takes into account 
the additional data for A — 1 shown in Fig. O Finally, 
there is a small structure at low magnetic fields in the 
A = 1 curves of Figs. HJc) and [SJa) which reflects the 
singular nature of the Heisenberg point (A = 1,H = 0). 

At the isotropic point of the antiferromagnetic XX Z 
chain, the quantity Th/H shows singular behavior for 
T -> if H < T. In fact, in the limit of H -> the ratio 
reduces to the temperature derivative of the zero-field 
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FIG. 2: (Color online) Cooling rate Vh for the Heisenberg 
chain (A = 1). Shown are exact diagonalization (ED) for 
TV = 20 sites, Quantum Monte Carlo (QMC) for N = 100 sites 
and the exact solution for the infinite system. The different 
panels are for T/J = 0.3025 (a), 0.1472 (b), and 0.1037 (c). 



FIG. 3: (Color online) Cooling rate Th for different values 
of the anisotropy. Shown is the exact solution for the infinite 
system. The different panels are for T/J = 0.05 (a), 0.5 (b), 
and 5 (c). 



susceptibility x(T) 



Urn ^ 



J_ fdx 

C H Ut 



(20) 



H 



which is known to show singular behavior. For A = 1 this 



function diverges like const, x [Tln(To/T)]~ 2 where To is 
a (non-universal) constant^ For A < 1 the function ([50]) 
diverges like const, x T 4K ~ 6 with Tomonaga-Luttinger 
parameter K = tt/(tt — 7) > l4i For A > 1 a divergent 
behavior like const. /T is observed. The strongest diver- 
gence is hence exhibited for A = 1. For finite magnetic 
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field H we see a non-divergent behavior for T — > 0. 

With increasing temperature, all features become 
broader. Note that the temperature T = 0.5 J shown 
in Fig. [3Jb) is higher than the value of H c \ at A = 2 
which explains why the sign change around H c \ in the 
A = 2 curve disappears for T = 0.5 J. For the even 
higher temperature T = 5 J shown in Fig.^c), the cool- 
ing rate Yh is positive for all H > (the zero at H = is 
enforced by the symmetry under H — > —H). Note that 
T = 5 J is bigger than _ff c2 itself for all cases shown in 
Fig. OJc) which explains why all features are washed out 
at this temperature. 

D. The case of infinitely large anisotropy: cooling 
rate for the Ising chain 

The case of infinitely large anisotropy (A — > oo at 
J i = J A = const.) corresponds to the Ising model with 



Hamiltonian 

N N 

H = jY,siSi+i- H J2 s *> ( 21 ) 

i=l i=l 

where the classical variables Sj take values ±1/2 and for 
simplicity reasons we drop the subscript of the coupling 
Jj in this section. It is well known that the Ising chain 
can be solved exactly and completely analyticall y by the 
transfer-matrix method (see for example Rcfs. 48.49) and 
even (dM/dT) H was investigated recently^ However, 
we are not aware of explicit results for the entropy S 
and normalized cooling rate T h of the Ising chain and 
therefore present them here. 

We start from the free energy per lattice site of the 
Ising chain: 



/ = ~ In (e- f3J l 4 (cosh(f3H/2) + ^sinh 2 (/Iff/2) + e /"^ 



(22) 



From this expression one can easily obtain simple analytic expressions for all thermodynamic functions of the system. 
In particular one obtains for the entropy per spin and the normalized cooling rate: 



d 2 f 



I ,. 1 ( '" \ - _ELUK 



-J/4(c + Q) + sH/2+ sH/2 + Je ' 3 
T(c + Q) 

l/2(ci?/2-.sJ/2)(c + Q) 2 



T V 9H J a 1/2 J 2 s 2 Q + 1/4 J 2 c(s 2 + Q 2 ) + H/2(cH/2 - Js)(c + Q) 2 ' 



where 



c = cosh(/3if/2) , s = sinh(/3#/2) , Q = \/smh 2 (f3H/2) + e$ J . 

r 



(23) 
(24) 

(25) 



The T = limit of the entropy (|23|) is generically S = 0, 

except for = ± J, where one finds S = In = 

0.4812 . . . (see also Ref . [5lh . This reflects the macroscopic 
ground-state degeneracy of the Ising model at the satu- 
ration field H c = JfiS Remarkably, the above transfer- 
matrix solution is closely related to the hard-dimer de- 
scription of certain highly frustrated one-dimensional 
quantum antiferromagnets i 12 ' 22 ! 52 ' 53 ' 54 

Fig. [4] shows the cooling rate dependence on the mag- 
netic field for the Ising chain obtained from the exact so- 
lution. The main difference from the XXZ case are the 
very sharp and pronounced positive and negative peaks 
at the critical value of the magnetic field. The magni- 
tude of these peaks grows rapidly with decreasing tem- 
perature. This behavior is a direct consequence of the 
anomalous zero-temperature entropy of the Ising chain 



at H = 

A bit further away from the QCP at H c = J, we can 
make contact with the argumentation of Ref. (Tj. The 
Griineisen ratio, which in our case is the cooling rate Th, 
shows divergent behavior close to the QCP at H c and 
changes sign when the magnetic field crosses it. For H ^ 
H c , the divergent behavior obeys the universal scaling 
lawl^ 



T H (T^0,H) = -G H 



1 



H-H r 



(26) 



where Gh is a universal amplitude. Detailed analysis of 
the Ising case yields the exact analytic form of the cool- 
ing rate at extremely low temperatures, r^ mg (T — > 0) = 
for T <g; \H — J\, i.e., the value Gh = —1 expected 



H-J 



for a Z2-symmetry in one dimension 
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FIG. 4: (Color online) Cooling rate Th of the Ising model 
(A — > oo limit of the XX Z model). Shown is the exact 
solution for the infinite system. The different panels show 
the cooling rate for a wider range of temperatures (a) and the 
behavior at low temperatures in the vicinity of the critical 
point (b). 



E. Entropy 

Finally, we take a look at the entropy S. As far as we 
are aware, there has been only one previous attempt of 
an exact computation 11 of the entropy and accordingly 
the isentropes of the XX Z chain in a magnetic field. 
Inspection of the critical fields H c \ and H C 2 indicates 
that the result of Ref. [ll| corresponds to an anisotropy 
A ps 4.75. The main technical difference between Ref. fill 
and the present work is that Ref. [ill worked with an 
infinite hierarchy of equations for A > 1 which need to be 
truncated while we have a closed finite set of equations 
for all values of A. In the case |A| < 1 the structure 
of the NLIE ([6]) is independent of A in contrast to the 
usual TBA equations, which allows to easily calculate the 
entropy as a function of the anisotropy. 

Fig. [5] shows the entropy and the isentropes for the 
spin-1/2 Heisenberg chain in the H-T plane. This result 
agrees with previous ED for N — 20 sitesJ^ However, 



0.1 



0.01 




FIG. 5: (Color online) Entropy per spin S(H, T) for A 
The isentropes are for S = 0.05, 0.1, . . . , 0.6. 



the ED results of Ref. [l2| suffered from finite-size effects, 
in particular for low temperatures and H < H C 2 — 2 J. 
By contrast, Fig. [5] shows results for the thermodynamic 
limit. Fig. shows a similar plot of the entropy and 
isentropes with varying anisotropy A but now at a fixed 
magnetic field H . The quantum phase transitions at H c i 
and H C 2 are reflected in Figs. [5] and [6] by minima of the 
isentropes as a function of H or A, or equivalently max- 
ima of the entropy at a low but constant temperature. 
The only exception is Fig. El^a) where one observes no 
such clear signature of the Kosterlitz-Thouless transition 
at A = 1 for H = 0. Before we discuss this case in more 
detail, it is useful to examine the low-temperature asymp- 
totics of the entropy S (T) at otherwise fixed parameters 
H and A. 

Between the two quantum critical points, that is for 
H c i < H < H c2 , the low-energy theory is a Luttingcr 
liquid (compare remarks at the beginning of Sec. |TTJ) . It 
is well known that the specific heat C of a Luttinger liquid 
is linear in TM Due to the relation C = T (d S/dT) (see, 
e.g., Ref. and because of S(T = 0) = 0, the entropy of 
a Luttinger liquid is identical to its specific heat and in 
particular also linear 



3r 



T for H cl < H < H, 



c2 ■ 



(27) 



where v is the velocity of the excitations. 

The cases H = H c i and H c2 ar e instances of a quantum 
phase transition in one dimension which preserves a C/(l)- 
symmetry. In this case, the universal low-temperature 
asymptotics is predicted to follow a square roo ^ 12 ' 16 i 17 



S oc \l — for H = H cl or H c2 . (28) 

Finally, in the gapped cases < H < H c \ or H > H C 2, 
we expect activated behavior 



S (x exp(-G/T) , 



(29) 
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(b) SforH/J=1.5 
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FIG. 6: (Color online) Entropy per spin S(A,T) for different 
values of H = (a), 1.5 J (b), and 2.5 J (c). The black 
lines are the isentropes for S = 0.05,0.1, . . . ,0.6. Panel (a) 
contains three additional isentropes with S = 0.005, 0.01, and 
0.02 which are shown by white lines. 



where G is the gap in the excitation spectrum. Closer 
inspection of the data shown in Figs. [5] and [6] indeed ver- 
ifies Eqs. P? |) , (|2g ]) . or ([25 ]) . respectively. Note that while 
the asymptotic behavior is S — > for T — > in all three 
cases, the decay is slowest exactly at the quantum phase 



FIG. 7: (Color online) Entropy per spin S of the XX Z chain 
at H = for A = 0.5, 1, and 1.5. For A = 0.5 and 1 we 
also show the Luttinger-liquid expression Eq. (|2T|) with the 
velocity v given by Eq. ((30} . 

transition, see Eq. (|28|) . This naturally yields a maxi- 
mum of the entropy S if the quantum phase transition 
is crossed by varying the parameters H or A at a fixed 
temperature T > 0. 

The case A = 1, H = is an exception to this gen- 
eral scenario. To first approximation, this point behaves 
like a Luttinger liquid. The fact that one is at a quan- 
tum critical point with marginally irrelevant operators 
gives rise to higher-order logarithmic corrections in the 
free energy and specific heati 29 i 30 Consequently, one also 
expects just higher-order logarithmic corrections to the 
entropy. To test this scenario, we can use the fact that at 
H = the velocity v which enters is known exactly 
(see, e.g., Refs. I39ll55[ ): 



with A = cos 7 as above. Fig. [7] shows that insertion 
of these values of v into the Luttinger-liquid expres- 
sion Eq. (|27[) yields indeed the correct low-temperature 
asymptotics of the entropy per spin at H = not only for 
A = 0.5 but also for A — 1. In fact, the higher-order log- 
arithmic corrections to the Luttinger liquid asymptotics 
P?)) which are expected at A = 1 and H — are so small 
that they have no visible effect. 

A gap G opens for A > 1 at H — 0, but because the 
phase transition is a Kosterlitz-Thouless transition ; 3 ^ 38 
this gap is exponentially small close to A = L 35 ' 36 Ac- 
cordingly, close to A = 1 one has to go to very low tem- 
peratures in order to observe the crossover from ([77]) to 
([29")). This is illustrated by the A = 1.5 curve in Fig. [7] 
In the concrete case A = 1.5, the value of the gap is 
G w 0.087 J. Accordingly, the exponential decay (|29]) 
can be observed only for temperatures T < 0.05 J while 
at higher temperatures the behavior of S(T) remains ap- 
proximately linear. 
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The combined effect of all these observations is that 
just a small kink develops in the low-temperature isen- 
tropes of Fig. \W[a) whose position shifts very slowly to 
the quantum critical point A = 1 for T — > 0. 

III. CONCLUSION AND PERSPECTIVES 

Motivated by recent measurements of the magnetic 
cooling rate in a spin-1/2 Heisenberg chain compound, 15 
we have presented an exact computation of the en- 
tropy and the magnetic cooling rate of the antiferromag- 
netic spin-1/2 XX Z chain in the thermodynamic limit 
N — > oo. Furthermore, we have performed complemen- 
tary numerical computations for the cooling rate of finite 
Heisenberg chains, namely exact diagonalization (ED) 
for small systems and QMC simulations for somewhat 
longer chains. We have demonstrated that we can obtain 
excellent approximations to the exact result with a com- 
bination of both numerical methods. On the one hand, 
this serves as a consistency check of our computations. 
On the other hand, we are now in a position to apply 
a combination of ED and QMC to some sign-free situa- 
tions, like the spin-1/2 Heisenberg model on the square 
and simple cubic lattices where the exact methods are no 
longer applicable. 

We have used the exact result for the entropy to illus- 
trate that field-induced quantum phase transitions give 
rise to maxima of the low-temperature entropy, or equiv- 
alently minima of the isentropes. This leads to cooling 
during adiabatic (de)magnetization processes where the 
lowest temperature is reached close to the quantum phase 
transition. As a consequence, we find a zero for the 
magnetic cooling rate at the phase transition and large 
positive (negative) values of the normalized cooling rate 
([2]) for magnetic fields slightly above (below) the critical 
field. 

The low-temperature asymptotics of the entropy S 



is exponentially activated in the gapped phases, linear 
in T in the gapless Luttinger liquid regions, and fol- 
lows the square-root behavior (|2"8|) at the field-induced 
quantum phase transitions. These asymptotic forms 
of S(T) are expected to be universal for field-induced 
phase transitions in one-dimensional systems with U(l)- 
sy mmetr y^ i 1 2 ' 1 6 1 1 7 but can be particularly clearly verified 
with the aid of an exact solution. 

The general features of the entropy should not depend 
on the specific choice of the magnetic field H as control 
parameter and indeed similar behavior is found as a func- 
tion of the exchange anisotropy A. An exception is just 
the quantum phase transition at H = 0, A = 1 with 
A as a control parameter. Because this is a Kosterlitz- 
Thoulcss transition, only weak signatures are observed in 
the finite-temperature entropy. 

Finally, closed-form expressions were derived in the 
Ising limit using the transfer matrix mcthod i 48 ' 49 We 
have observed remarkably large magnetic cooling rates 
close to the field-induced critical point of the Ising 
chain. In fact, the transfer-matrix solution is closely 
related to a low-energy description of highly frustrated 
one-dimensional quantum ant if err omagnet s ) 1 2 ' 22 1 5 2 ' 53 1 54 
where enhanced cooling rates are observed as well. 
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